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In radio bands, the study of compact radio sources can be affected by propagation effects through 
the Interstellar Medium. These are usually attributed to the presence of turbulent intervening 
plasma along the line of sight. In this talk, two of such effects are presented. The line of sight 
of B 2005+403 passes through the heavily scattered region of Cygnus, which causes substantial 
angular broadening of the source images obtained at frequencies between 0.6 GHz and 8 GHz. 
At higher frequencies, however the intrinsic source structure shines through. Therefore, multi- 
frequency VLBI observations allow to study the characteristics of the intervening material, the 
source morphology and the "coupling" of them in forming the observed image. 
This article is based upon the published paper of Gabanyi et al. For more details see m\. 
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1. Introduction 

The propagation of radio waves through the ionized interstellar medium causes several effects, 
such as Faraday rotation and depolarization of polarized emission, dispersion of pulsar signals, 
scatter broadening of compact radio sources and intensity fluctuations caused by diffractive (DISS) 
and refractive (RISS) interstellar scintillation (ISS). The detailed study of scatter broadening in 
compact radio sources can lead to a better understanding of the interstellar medium. 

The other prominent propagation effect addressed here is the so-called Intraday Variability 
(IDV, and the question of how it is related to the interstellar scintillation. IDV surveys show 
that a large fraction (up to 30 %) of all compact flat spectrum radio sources show this effect, which 
is characterized by variability amplitudes of up to 20 — 30 % and variability timescales ranging 



from less than one hour to several days (e.g. [10], [12] and references therein). If interpreted 



via source intrinsic incoherent emission processes, such short variability timescales imply - via 
the light travel time argument - apparent brightness temperatures of Tg = 10 16 K — 10 19 K, far in 
excess of the inverse-Compton limit ([§]). With the assumption of relativistic Doppler-boosting the 
brightness temperatures can be reduced, which however would imply uncomfortably large Doppler 
boosting factors of 8 ~ 20 — 200. 

Alternatively, IDV in the radio-bands is also explained extrinsically via scintillation of radio 
waves in the turbulent interstellar medium (ISM) of our Galaxy (e.g. [|l^]). The main problem in 
this interpretation is that it cannot explain the observed radio-optical broad-band correlations seen 



in at least some sources ([|14|]). The recent detection of diffractive ISS in the source J 1819+3845 
leads to micro-arcsecond source sizes and brightness temperatures of ~ 10 14 K, which again require 
Doppler boosting factors of 5 ~ 100 ([|j"3|]). It is therefore likely that the IDV phenomenon involves 
both, a combination of source intrinsic and extrinsic effects (e.g. [[fl|]). 

We present measurements for the scatter broadened quasar B 2005+403, combining Very Long 
Baseline Interferometry (VLBI) and flux density variability measurements obtained at various fre- 
quencies. The scatter broadening observed in the VLBI images at the lower frequencies is measured 
and yields further constraints to the properties of the intervening ISM. Additional parameters for 
the ISM are obtained from the IDV monitoring of the source, performed with the Effelsberg 100 m 
telescope. 

B 2005+403 is a flat spectrum quasar (G!o.3/5GHz = 0-3, with S ~ v a ; [jj) at a redshift of 
z = 1.736 ([§]). It is located close to the Galactic plane at 1 = 76.82°, 6 = 4.29° behind the Cygnus 
super-bubble region. Earlier studies showed that interstellar scattering affects the VLBI image 
of the source, causing angular broadening at frequencies below 5 GHz (see [Q] and references 
therein). The high value of the scattering measure (SM = / L Cft(s)rfs ! ), derived for the line of 
sight of B 2005+403 by reflects the strong influence of the interstellar medium. 

2. VLBI data and data reduction 

The VLBI observations of B 2005+403 were performed at frequencies between 1.6 GHz and 
43 GHz. 



'The Scattering Measure is the path integral over the coefficient Cj| of the electron density fluctuation wavenumber 
spectrum. 
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Figure 1: The measured angular size plotted 
versus observing frequency. The solid black line 
represents a power law fit to the data in the range 
from 0.6 GHz to 8 GHz. The dashed blue line 
represents a power law fitted to the data in the 
range from 15 GHz to 43 GHz. Black circles de- 
note data from the literature ([§], [Q]), red trian- 
gles denote our data. 
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Table 1: The measured angular sizes of 
B 2005+403 at different frequencies. 



After correlation at the VLBI correlators in Socorro (NRAO) or Bonn (MPIfR), a priori am- 
plitude calibration using system temperature measurement and fringe-fitting were performed using 
the standard AIPS analysis tasks. Editing, phase and amplitude self-calibration and imaging were 
performed using AIPS and the Caltech DIFMAP packages. 

To study the scattering effects following [||], we measured the total angular source size as a 
function of observing frequency by fitting the visibility data of B 2005+403 at all frequencies with 
one Gaussian component within DIFMAP. 

Additional size measurements of B 2005+403 were used from []|] and g]. 

3. Propagation effect I: Scatter broadening of B 2005+403 

The angular source size of B 2005+403 versus frequency is plotted in Fig. [j], using our 
data and including the previous measurements from the literature. The corresponding values 
are given in Table [jj To quantify the scatter broadening, which dominates at lower frequencies, 
a power law was fitted to the size-frequency relation in the frequency range from 0.67 GHz to 
8 GHz, yielding d = (77.1 ±4.0) • (v/1 GHz)-( L91±0 05 ) mas. This fit is shown as the black line 
in Figure [j]. (Excluding the data point at 8 GHz does not change the slope significantly, yielding: 
d = (77.7 ±4.0) • (v/1 GHz) _1 - 92±0 06 .) Above 8 GHz, the extrapolated scattering size becomes 
smaller than the measured source size. This indicates that towards higher frequencies scattering 
effects are less dominant and that the intrinsic structure of the sources shines through. The differ- 
ences between the extrapolated scattering size and the measured source size at all three frequencies 
(15 GHz, 22 GHz, and 43 GHz) are of the order of 0.4 mas. 

To characterize the frequency dependence of the intrinsic source size, a power law was fitted to 
the size-frequency relation above 8 GHz, yielding: int = (1.7 ±0.4) • (v/1 GHz)~(°- 31±007 ) mas. 
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Source name 


(s) [Jy] 


a[Jy] 


m[%] 


Y[%] 


x r 2 


V = 1.6 GHz, mo = 0.25 % 


B 2005+403 


2.430 


0.024 


1.01 


2.93 


5.593 


NGC 7027 


1.906 


0.004 


0.19 





0.197 


B 2021+614 


2.193 


0.006 


0.25 





0.365 


V = 5 GHz, m 


= 0.20% 










B 2005+403 


2.905 


0.013 


0.45 


1.18 


3.107 


NGC 7027 


5.489 


0.011 


0.20 





0.621 



Table 2: The variability parameters of B 2005+403 and the secondary calibrators at 1 .6 GHz (top) and 5 GHz 
(bottom). Col. 2 shows the mean flux density in Jy, col. 3 its standard deviation, col. 4 the modulation index, 
col. 5 the noise-bias corrected variability amplitude and col. 6 the reduced Xi- For a detailed definition of 
these values see [ p~0|| . 



The dashed blue line in Fig. [j] shows this fit. With this line, it was possible to obtain an upper 
limit on the intrinsic source size at lower frequencies where direct size measurements with VLBI 
are not possible. Hence a lower limit on the brightness temperature could also be calculated. 
Below 8 GHz, the intrinsic source sizes are in the range of 1 — 1.5 mas. The corresponding lower 
limits to the brightness temperatures at 1.6 GHz, 5 GHz and 8 GHz are T B > (0.6 ±0.5) • 10 12 K, 
T B > (0.14 ±0.09) • 10 12 K, and T B > (0.08 ±0.04) • 10 12 K, respectively. These numbers are in 
accordance with typical brightness temperatures measured with VLBI and neither strongly violate 
the inverse-Compton limit nor do they indicate excessive Doppler-boosting. 

The size of the scattering disk at a given frequency can be estimated via the deconvolution 
formula: scat = J 2 bs — d 2 nV From the power law fit, the scattering size at 1 GHz Qy G Hz is (77. 1 ± 
4.0) mas. Following [16] and assuming a Kolmogorov turbulence, the scattering measure can be 
calculated. The derived SM of (0.43 ± 0.04)m~ 20 / 3 kpc is in good agreement and consistent with 
the previous measurement of [S], though with much improved accuracy. 



3.1 Propagation effect II: IDV behaviour of B 2005+403 

In November and December 2003 the flux-density variability of B 2005+403 was monitored 
with the Effelsberg 100 meter radio telescope with high time resolution. In Table Q the results 
from these measurements are summarized: To characterize the variability amplitudes and their 



significance, we follow the methods described in Q10|]. Column 2 of the table gives the average flux 



density, col. 3 the rms standard deviation, col. 4 the modulation index (m = 100 • a/{S)), col. 5 the 



variability amplitude (defined as Y = 3y m 2 — m^, where mo is the modulation index of a source 
regarded to be stationary during the observation) and col. 6 the reduced % 2 . A value of Y is only 
given, if the # 2 — test gives a higher than 99.9 probability for significant variations. 

The variability amplitude of B 2005+403 decreases with frequency. Formally we measure a 
modulation index of m = 1.01 % at 1.67 GHz and of m = 0.45 % at 5 GHz. At this frequency the 
detection of IDV is marginal. We note that the model of refractive interstellar scintillation in the 



weak regime predicts a decrease of the modulation index with increasing frequency ([15]). This is 
consistent with our observations. 
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The source size and the apparent brightness temperature can be determined from the interstel- 
lar scintillation model, assuming that the main reason for the observed variability is the motion 

of the Earth through the scintillation pattern. In this scenario the variability time-scale in days is 

-l 



(@]):^ = 18-0 scat -L kpc . V5Okms _ 

Cordes and Lazio modeled the electron density turbulences of the Milky Way (NE2001, [J3j] ) - 
They assumed a distance of 2.35 kpc to the scattering "clump" responsible for the angular broaden- 
ing of B 2005+403. The characteristic variability time-scales derived from the light curves is less 
than a day in B 2005+403. A screen located at kpc distance cannot explain the observed short vari- 
ability time-scale. With a scattering size of scat = 30 mas at 1.6 GHz, a screen distance of 2.35 kpc 
and typical Galactic velocities of < 220km/s, one would expect to see variations on time-scales 
> 288 days. 

To reproduce a variability time-scale of < 1 day at 1.6 GHz within the aforementioned model, 
the following constraints for the scattering size and the screen distance are obtained: e ffL < 
2.2 • 10~ 2 mas kpc. For the measured scattering size of ~ 30 mas, this leads to an unreasonably 
nearby screen of < 1 pc. The only way out of this dilemma is a smaller scattering size. Adopt- 
ing a minimum screen distance of at least 10 pc as required for the interpretation of the ultra fast 
scintillators ([13] and references therein), one obtains scat < 2.2 mas. 

An upper limit on the distance of the screen is obtained from the restriction of the source size 
via the inverse-Compton limit of Tg C = 10 K. The requirement that the brightness temperature 



be lower than this leads to a source size of 8[ nt = J 1.77 • \0 12 S v v^ 2 T Ji l 8^ l (l +z) > 0.6mas, 
where Si .6GHz = 2.4 Jy is adopted from Table ^| and 8 ~ 10 for the Doppler factor and V is in 
the units of GHz. Approximating scat with this size and a relative screen velocity in the range 
(50 to 220) km/s, leads to an upper limit of the screen distance in the range (9 to 41) pc. 



4. Conclusion 



The observed IDV in B 2005+403 can be explained with a screen located at a distance less than 
or equal to 41 pc that can be characterized by a scattering size of 0.6 mas < scat < 2.2 mas. The 
corresponding scattering measure at 1.6 GHz then is in the order of 10~ 3 m~ 20 / 3 kpc. These values 
are considerably lower than the scattering measure of the more distant screen, which is thought to 
be responsible for the scatter broadening in B 2005+403. Therefore the observed IDV is not caused 
by this distant screen. Most likely, one has to assume multiple scattering by at least two spatially 
and physically very different plasma screens. In this scenario, the first screen leads to a significant 
scatter broadening of the source image, which then by the second screen leads to only very weak 
scattering due to large quenching effects (Source >> scat ). This quenched scattering ([|l^]) by the 
second screen can explain the relatively low variability amplitudes of < 1 %, observed in the IDV 
experiments at 5 GHz and 1.6 GHz. Quantitatively, this can be verified using equation (20) of [pi], 
which relates the variability index, the scattering measure, the effective source size and the screen 
distance. With the parameters from above, a modulation index of m < 1 % is obtained, in good 
agreement with the observations. 
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